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I _ Abstract: We calculate the complete next-to-leading order (NLO) QCD corrections to 

the Z^H^ production in association with a jet at the LHC. We study the impacts of the 
NLO QCD radiative corrections to the integrated and differential cross sections and the 
dependence of the cross section on the factorization/renormalization scale. We present 
the transverse momentum distributions of the final Z^-, Higgs-boson and leading-jet. We 
find that the NLO QCD corrections significantly modify the physical observables, and 
obviously reduce the scale uncertainty of the LO cross section. The QCD K-factors can 
be 1.183 and 1.180 at the y/s = 14 TeV and ^/s = 7 TeV LHC respectively, when we 
adopt the inclusive event selection scheme with p^j = 50 GeV, mn = 120 GeV and 
fi = fj,r = fJ'f = P'0 = \{'^z + i^h)- Furthermore, we make the comparison between the 
two scale choices, = /xq and ^ = fii = \{E^ + + J2j Et^)i and find the scale choice 
fi = fii seems to be more appropriate than the fixed scale fj, = fj,Q. 

Keywords: Higgs-boson production, NLO QCD Corrections, Large Hadron Collider 
(LHC) 
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1. Introduction 

The standard model (SM) is a non-Abelian gauge theory and most of its predictions have 
been perfectly confirmed by the present precise experimental data except the Higgs-boson, 
which explains the mechanism of mass generation and is believed to be responsible for 
the breaking of the electroweak symmetry [1-4]. Therefore, one of the primary tasks for 
the CERN Large Hadron Collider (LHC) is to discover the Higgs-boson. Correspondingly, 
major efforts have been concentrated on devising method for Higgs-boson search. The four 
LEP collaborations have established the lower bound of the SM Higgs mass as 114.4 GeV 
at the 95% confidence level [5]. The CDF and DO experiments at the Fermilab Tevatron 
have ruled out the SM Higgs boson with mass between 156 and 177 GeV at 95% confidence 
level (CL) [6]. The ATLAS and CMS experiments at the LHC have excluded most of the 
niH ranges 146 — 466 GeV and 145 — 400 GeV in latest reports [7] and [8] respectively, at 
95% CL. The present LHC data are not yet sensitive to the region around itih ~ 120 GeV 
favored by the SM fit to electroweak precision data [9]. Not long ago, the ATLAS and CMS 
experiments at the LHC have provided the upper limit of the SM Higgs mass as 130 GeV 
and 127 GeV at 95% CL respectively, and there are several Higgs like events around 
the locations of mn ~ 126 Gey(ATLAS) and mn ~ 124 GeV (CMS) [10] [11]. Further 
searching for Higgs boson and studying the phenomenology concerning its properties are 
still the important tasks for the present and upcoming high energy colliders. 

There are several Higgs search channels, such as: gluon fusion {gg — > H^) induced 
mainly by a top-quark loop [12, 13], weak boson fusion (WBF) [14], top-quark associ- 
ated production (ttH^) [15], and weak boson associated production {VH°{V = W^,Z^)) 
[16-18]. Among these channels, the VH'^ {V = W^,Z^) associated production processes 
are promising channels at low Higgs mass region. Although the production rates of the 
VH^ {V = W^,Z^) associated production processes are several times lower than that of 
the gluon fusion channel, these association productions become available for Higgs search- 
ing by using the techniques of jet reconstruction, b-tagging and lepton identification [19- 
22]. Furthermore, the VH^ {V = W^,Z^) associated production processes also provide 
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unique information on the couplings between Higgs and vector gauge-bosons. Recently, 
the calculations of the QCD 0{as) and electroweak 0{aew) corrections to the Higgs pro- 
duction processes pp/pp W^H^/Z^H^ + X aX the Tevatron and LHC were presented 
in Refs. [23-27] and Ref.[18], respectively. The NNLO QCD corrections to the SM Higgs- 
boson production processes associated with a vector boson at hadron colliders have been 
calculated in Ref.[16, 28]. In Ref.[28] it concluded that the NNLO QCD corrections to the 
Z^H^/W^H^ production processes at the Tevatron and the LHC increase the cross sec- 
tions by the order of 5 — 10% , while the electroweak 0{aew) corrections have been turned 
out to be negative and about —5% or —10% depending on whether the weak couplings 
are derived from or aewitn^z) ^respectively [18]. Therefore, after the inclusion of these 
corrections, the remaining uncertainties for the Z^H^ /W^H^ production processes should 
be dominated by factorization/renormalization scale dependence and parton distribution 
functions. 

As we known, the experimental environment at hadronic collider is extremely compli- 
cated. The VH^ {y = W^,Z^) associated production signals at the LHC are normally 
accompanied by multi-jets in final state. Therefore, it is very necessary to have a good 
understanding of these multi-body final state processes. That requires sufficiently precise 
predictions for the VH^ {V = W^, Z^) associated production signals and their backgrounds 
with multi-jets in final state which cannot entirely be separated in experimental data. Actu- 
ally, the inclusive VH^ (V = W^, Z^) associated production signals include any number of 
additional jets unless otherwise stated. In this sense the VH^ + jet {V = W^, Z^) produc- 
tions are part of the inclusive VH^ {V = W^, Z^) productions, and theoretically VH^ +iet 
production at the next-to-leading order (NLO) QCD is part of the VH^ production pro- 
cess at the QCD next-to-next-to-leading order (NNLO). The process pp — t- W^H'^j + X at 
hadron colliders including the complete NLO QCD corrections, which are part of the NNLO 
QCD correction to the pp/pp — )• W^H^ + X process, have been calculated in Ref.[29]. 

In this paper, we present the calculations for the process pp — )• Z^H^j + X at the 
LHC up to the NLO in the QCD. The paper is organized as follows: In Sec. 2, we describe 
the calculations of the tree-level cross section and the QCD NLO corrections for the pp — )• 
Z^H^j + X process at the LHC. The numerical results and discussions are presented in 
Sec. 3. A short summary is given in Sec. 4. Finally, we provide some analytical expressions 
for the amplitudes of the LO and real emission partonic processes in Appendix. 

2. Calculations 
A. LO calculations 

The LO cross section for the pp — t- Z^H^j+X parent process involve following partonic 
processes: 

q (pi) + q {P2) ^ Z^ (p3) + (pa) + 9 fe), (2.1) 

q {Pi) + 9 {P2) ^ ^° (P3) + {pa) + q (ps), (2.2) 

q (Pi) + 9 {P2) ^ Z^ (ps) + (pa) + q{p5), (2.3) 
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(4) (5) (6) 



Figure 1. The generic LO Feynman diagrams for the partonic processes qq — >■ Z^H^g, qg — >■ Z^H^q 
and qg — >■ Z^H^q. (1) and (2) are the LO diagrams for the partonic process qq — )• Z^H^g, (3) and 
(4) for the qg Z^H^q, (5) and (6) for the qg Z°H°q, where q = u,d, c, s, b. 



where q = u,d,c,s,b and pi {i = 1,...,5) represent the four-momenta of the incoming 
partons and the outgoing Z^, H^, jet, respectively. There are 6 tree- level Feynman dia- 
grams for the partonic processes of the Z^H^+iet production, (2.1)-(2.3), shown in Fig.l. 
There Figs. 1(1-2) are the LO diagrams for the partonic process qq — t- Z^H^g, Figs. 1(3-4) 
for qg — t- Z^H^q and Figs. 1(5-6) for qg — ?• Z^H^q. They can be grouped into two different 
topologies. Figs. 1(1, 2, 3, 5) belong to u(t)-channel, while Figs. 1(4, 6) to s-channel. Actually, 
the partonic processes of the Z'^i^'^+jet production at the LHC are related to the amplitude 
of — )• Z^H^qqg by crossing symmetry at the LO. 

The LO cross section expressions for the partonic processes qq — t- Z^H^g, qg — t- Z^H^q 
and qg — >• Z^H^q have the forms respectively as: 

spin spin 

(2.4) 

where A4|^q, -^io -^io amplitudes of the corresponding tree-level diagrams 

for the partonic processes qq — t- Z^H^g, qg — t- Z^H^q and qg — t- Z^H^q (q = u, d, c, s, b) 
shown in Fig.l. The factors j and ^(^) in Eqs.(2.4) are due to spin- and color- averaging 
for the initial partons, respectively, s is the partonic center-of-mass energy squared. The 
summations in Eqs.(2.4) are taken over the spins and colors of all the relevant particles in 
these partonic processes, dfig' {kl = qq,qg,qg) in Eqs.(2.4) is the three-body phase space 
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element expressed as 



5 5 



1=3 j=3 



The LO total cross section for the parent process pp — t- Z^H^j + X can be expressed 



as 

aLo{AB{pp)^Z^H''j + X) 

= ^ j dxAdxB Gk/A{xA,f^f)Gi/B{xB,fJ'f)^LoixAXB,fJ'f) + iA^ B) . (2.6) 
kl 

Here kl = qq,qg,qg, {q = u,d,c,s,b), AB = pp and is the factorization scale. xa{xb) 
describes the probability to find a parton k{l) in proton A{B) defined as 

Pi P2 „^ 

XA = ^, XB = ^. (2.7) 

Here Pa and Pb are the four-momenta of the corresponding protons. Gi.^iyA(B) ^'^^ the 
LO parton distribution functions (PDFs) for parton k[l) in a proton. 
B. NLO QCD corrections 

Throughout our calculations, we take the 't Hooft-Feynman gauge except when we 
verify the gauge invariance. In the NLO QCD calculations we adopt the dimensional reg- 
ularization (DR) method in D = 4 — 2e dimensions to isolated the UV and IR singularities 
and the modified minimal subtraction (MS) scheme to renormalize the colored fields. The 
one-loop diagrams are essentially obtained from the tree-level diagrams of related partonic 
processes (2.1)-(2.3) and generated by means of the FeynArts3.5 package [30]. The am- 
plitudes are further analytically simplified by the modified FormCalc programs [31]. The 
reduction of a tensor integral to the lower-rank tensor and further to scalar integral is 
done with the help of the Loop Tools library [31, 32] and the FF package [33]. There the 
dimensionally regularized 3- and 4-point integrals with soft or collinear singularity have 
been added to this library [34, 35]. The output amplitudes are numerically evaluated by 
using our developed Fortran subroutines for calculating N-point integrals and extracting 
the remaining finite e x ^ terms. 

The NLO QCD corrections to the parent process pp — ?• Z^H^j + X involve following 
components: 

• The virtual contribution to the partonic channels (2.1), (2.2) and (2.3) from the QCD 
one-loop and the corresponding counterterm diagrams. 

• The contribution of the real gluon emission partonic processes. 

• The contribution of the real light- (anti) quark emission partonic processes. 

• The corresponding contribution of the PDF counterterms. 
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The one-loop QCD contribution to the pp — t- Z^H^j + X process contains both ultra- 
violet (UV) and soft/collinear infrared (IR) singularities, but the total virtual correction 
contributed by both the one-loop QCD and the corresponding counterterm diagrams is 
UV finite after performing the renormalization procedure. Nevertheless, it still contains 
soft/collinear IR singularity which can be canceled by adding the contributions of the real 
emission partonic processes and the corresponding PDF counterterms. 

The virtual corrections to the partonic processes qq — )• Z^H^g, qg — t- Z^H^q and 
qg Z^H^ q can be expressed as 



4 9 2s 

spin 



dnf, 



color 



= 1^ ^ 27^e Mlf'^^Mf''^ dnf''\ (2.8) 



spin 



where A^^q, -^^q a-nd A4^q are the LO matrix elements of the partonic processes qq — )• 
Z^H^g, qg Z^H^q and qg Z^H^q separately, and Mf, M^^ and M^^ are the NLO 
QCD virtual matrix elements of the corresponding partonic processes. 

The real gluon/light-(anti) quark emission partonic processes are obtained from — ?• 
Z^H^ggqq and — t- Z^H^qqq'q' by all possible crossings of (anti)quarks and gluons into the 
initial state. All the real gluon/light-(anti)quark emission partonic processes are presented 
in the following: 



(1) a{pi) 9{P2) - 


^°(P3) 


h^{pa) q{pb) q{pfs), 


(2.9) 


(2) q{pi) q{p2) - 


^ Z\pz) 


h^{pa) g{pb) g{pQ), 


(2.10) 


(3) q{pi) g{p2) - 


^ Z'ips) 


H^{pa) q{P5) giPe), 


(2.11) 


(4) q{pi) g{p2) - 


^ Z\ps 


H^{pi) q{P5) giPd), 


(2.12) 


(5) q{pi) q{p2) - 


■> Z^iPs) 


H\pa) q'{P5) q'iPe), 


(2.13) 


(6) q{pi) q'{p2) - 


Z\pz 


) H^iPi) q{P5) q'iPe), 


(2.14) 


(7) q{pi) q'{p2) - 


Z\p^ 


) H^{P4) q{P5) q'iPe), 


(2.15) 


(8) qipi) q'{p2) - 


Z\p^ 


) H^{P4) q{P5) q'{P6), 


(2.16) 



where q, q' = u, d, c, s, b. There are totally 95 real emission partonic channels, and these 
tree-level partonic processes contain both soft and collinear IR singularities. After the 
summation of the renormalized virtual corrections with all the real partonic emission cor- 
rections, the result is soft IR-safe, but still contains remained collinear divergence. It will 
be totally IR-safe when we include the contributions from the collinear counterterms of the 
PDFs. 

As a demonstration we show in Figs.2(a-h) the tree-level Feynman diagrams for the 
real emission partonic process qg — t- Z^H^qg. We adopt the two cutoff phase space slicing 
(TCPSS) method [36] to isolate the IR singularities by introducing two cutoff parameters 
6s and 6c- An arbitrary small 6s separates the four-body final state phase space into 
two regions: the soft region {Eq < 6sVS/2) and the hard region (£^6 > 6sV^/2). The 
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Sc separates hard region into the hard cohinear (HC) region and hard noncohinear (HC) 
region. The criterion for separating the HC region is described as foUows: the region for real 
gluon/hght-(anti) quark emission with siq (or S15, §25, §26, sse) < (where Sij = {pi+PjY) 
is called the HC region. Otherwise it is called the HC region. Then the cross section for 
each of the real emission partonic processes (2.9) can be written as 

a =a +a =a +a +a . [^■^'j 



After combining the renormalized virtual corrections with the contributions of the real 
gluon/light-quark emission processes and the PDF counterterms 5Cg(^gyp together, the UV 
and IR singularities are exactly vanished. These cancelations can be verified numerically 
in our numerical calculations. 

3. Numerical results and discussions 

In this section we present and discuss the numerical results for the LO and QCD NLO 
corrected observables for the pp — t- Z^H^j + X process at the early {^/s = 7 TeV) and 
future {^/s = 14 TeV) LHC. We take the CTEQ6L1 PDFs with a one-loop running 
in the LO calculations and the CTEQ6M PDFs with a two-loop running in the NLO 
calculations, separately [4, 37]. The number of the active flavors is taken as Nj = 5 and 
the QCD parameters are set as = 165 MeV and A^^^ = 226 MeV for the LO and 
NLO calculations, respectively. For simplicity we deflne the factorization scale and the 
renormalization scale being equal (i.e.,/i = fif = fir) and take = fiQ = ^{rriH + mz) by 
default unless otherwise stated. The CKM matrix is set as a unit matrix and the weak 
mixing angle is obtained by = m^^/m^. We neglect up-, down-, charm-, strange- 
and bottom-quark masses (i.e., m„ = = rric = nig = mi, = GeV). Throughout our 



calculations, we take mu = 120 GeV by default and the other input parameters are chosen 
in accordance with [38], 

a(m|)^i = 127.916 GeV, mw = 80.399 GeV, mz = 91.1876 GeV, mt = 172.0 GeV. 

(3.1) 

As we know that the final state of the pp — )■ Z^H^j + X process contains only one 
jet at the LO, while it involves both one-jet and two-jet events in final state up to the 
QCD NLO. The NLO virtual correction, the real soft emission and hard collinear emission 
processes include one-jet events, while the real hard noncollinear gluon/light-(anti) quark 
emission processes include two-jets events. In our calculations we apply the jet recombina- 
tion procedure of Ref.[39] in the definition of the tagged hard jet with R = 1. That means 
when two proto-jets in the final state satisfy the constraint of -sj Ay^ + /S.(jP' < R = 1, 
where Ay and are the differences of pseudo-rapidity and azimuthal angle between the 
two proto-jets, we merge them into one new 'jet' and call it as one-jet event, the new 
'jet' four- momentum is defined as Pij,^ = Pi,fj, + Pi,fi- After applying the jet recombination 
procedure to the proto-jet events of the pp — > Z^H^j + X up to the QCD NLO, we obtain 
the one-jet events and two-jet events. As shown in Ref.[40], the QCD NLO corrections 
to the total cross sections at hadron colliders are generally very large and could destroy 
the convergence of perturbative description, due to the fact that the contributions from 
NLO real emission subprocesses are taken into account. In order to know how to get the 
modest corrections and reduce the scale uncertainty, we adopt both the inclusive and ex- 
clusive two-jet event selection schemes for comparison. In the following we present detailed 
cr iter ions for selecting the one-jet and two-jet events. 

(1) For the one-jet events, we collect the events with the constraint on the jet as 
> pft = 50 GeV. 

(2) For the two-jet events, we treat the two jets either inclusively or exclusively. The 
inclusive and exclusive two-jet event selection schemes are declared as follows: 

• In the inclusive two-jet event selection scheme (Scheme I) both one- and two-jet events 
are included, and the constraint of piji^ > p^j = 50 GeV applied on the leading-jet 
but not on the second-jet, where the leading-jet and the second-jet are characterized 

1 leading— jet ^ second— jet 
by > prp . 

• In the exclusive two-jet event selection scheme (Scheme II), the one-jet events with 
p^^ > p^j = 50 GeV are accepted, while the two-jet events with p^^'^°^^ '^^^ > p^j = 
50 GeV are rejected [41-44] . 

The verifications of the correctness of our calculations are made in the following ways: 

(1) The LO calculations for the process pp ^ ug ^ Z^H^u + X with ^/s = 14 TeV are 
performed by using the FeynArts3.5/FormCalc6.0 packages and CompIIEP-4.4p3 program 
[45], and applying the Feynman and unitary gauges, separately. All the results are in good 
agreement within the statistic errors. 

(2) The virtual correction and the real gluon/light-(anti)quark emission correction to 
the pp —7- Z^H^j + X process at the LHC were evaluated twice independently based on 
different codes, and they yield results in mutual agreement. 
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(3) The UV and IR safeties are verified numerically after combining all the NLO QCD 
contributions. 

(4) The independence of the NLO QCD correction to the pp — )• Z^H^j + X process 
on the soft cutoff 6s and collinear cutoff 6c is proofed. There we apply the Scheme I for 
two-jet event selection, and take ^/s = 14 TeV, n = fiQ, ijih = 120 GeV and 6c = 6s/50. 
We find the total NLO QCD correction AaNLO does not depend on the arbitrarily chosen 
value of 6s and 6c within the calculation errors. In the following numerical calculations, we 
Rx6s = 5x 10-^ and 4 = 6s/50. 

In Figs.3(a,b) we depict the dependence of the LO, NLO QCD corrected cross sections 
and the corresponding K-factor for the pp — ?• Z^H^j + X process on the renormaliza- 
tion/f actor ization scale (n) at the ^/s = 14 TeV and -y/i = 7 TeV LHC, by adopting the 
inclusive and exclusive event selection schemes, separately. Figs. 3(a) and (b) show that the 
LO curves go down quickly with the scale running from 0.2 /xq to 5fXQ, while the renormal- 
ization/factorization scale dependence is obviously reduced by the NLO QCD corrections 
in both inclusive and exclusive event selection schemes. The scale uncertainty is defined 
^ V = . where /. e [0.2^0, S/xol. we can find that at the = 14 TeV 

LHC with the inclusive event selection scheme (Scheme I), the scale uncertainty rj is re- 
duced from 25.9% (LO) to 11.7% (NLO) as shown in Fig.3(a). At the = 7 TeV LHC 
with the inclusive event selection scheme, the scale uncertainty is reduced from 35.8% 
(LO) to 10.7% (NLO) as shown in Fig. 3(b). Figs.3(a,b) show also when we apply the ex- 
clusive event selection scheme (Scheme II), the NLO QCD corrections can reduce the scale 
uncertainty from 35.8% (LO) to 5.9% and 9.7% (NLO) with ^/s = U TeV and 7 TeV, 
respectively. We can also see that the K-factor varies between 0.996 (0.69) and 1.46 (1.55) 
at y/s = 14 TeV (7 TeV) by adopting the inclusive event selection scheme, and between 
0.50 (0.35) and 1.27 (1.42) at ^ = U TeV (7 TeV) by adopting the exclusive event se- 
lection scheme in the plotted range of fi. Comparing the curves in Figs.3(a,b), we can see 
the different line shapes for the two curves labeled with 'NLO (I)', which are for the cross 
sections obtained by adopting the inclusive scheme (Scheme I) at the ^/s = 14 TeV and 
-y/s = 7 TeV LHC, respectively. That difference comes from that the NLO QCD correction 
component {<J^^) of the non-collinear real hard emission correction at the = 14 TeV 
LHC is proportional to ^^(/u), and its contribution to integrated cross section is much larger 
than that at the y/s = 7 TeV LHC. Therefore, the line shape for t/s = 14 TeV in Fig. 3 (a) 
shows the cross section adopting the inclusive scheme is going down with the increment of 
the scale. 

In Table 1, we list the numerical results for the LO, NLO QCD corrected cross sections 
and the corresponding K-factor (K = ) for the process pp — t- Z^H^j + X at the LHC 
by applying the Scheme I and Scheme II with = /^o, ^ = /^i = ^{E^ + + Eip) 
and ^/s = 14 TeV, 7 TeV , separately. From this table, we can see that the LO (ctlo) and 
NLO QCD corrected cross sections by adopting the inclusive and exclusive event selection 
schemes (fj^^^Q, '^^{^q) are all sensitive to the scale choices, and the difference between the 
NLO QCD corrected cross sections obtained by taking fi = fiQ and fi = ni respectively, is 
about 3%, which is much smaller than the discrepancy of LO cross sections with these two 
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Figure 3. The LO, NLO QCD corrected cross sections and the corresponding K-factor {K{^) = 
'^NLo{^J■)/(^Loi^J')) versus the factorization/renormalization scale /i for the pp — > Z^H'^j + X pro- 
cess at the LHC by adopting the inclusive scheme (Scheme I) and exclusive scheme (Scheme II), 
separately, (a) at the Vs = 14 TeV LHC. (b) at the = 7 TeV LHC. 



pp Z^H^i + X 


fi{GeV) 


<TLo{fb) 






Kin 




^ = 14 TeV 


pto 


196.25(1) 


232.2(2) 


178.7(2) 


1.1832 


0.9106 


/^i 


175.03(1) 


226.4(6) 


184.9(6) 


1.2933 


1.0564 


^ = 7 TeV 




53.116(3) 


62.67(7) 


52.50(7) 


1.1799 


0.9884 


/^i 


46.022(3) 


60.86(5) 


53.71(5) 


1.3225 


1.1670 



Table 1. The numerical results for the LO, NLO QCD corrected cross sections and the 
corresponding K-factor {K = ^^) with p^!^ = 50 GeV, tuh = 120 GeV for the pp 
ZOffOj X process at the LHC by taking /i = = ^{mn + mz) = 105.594 GeV and 
H = Hi = \{E^ + E^ + Sjet ^T^*) E^nd applying the Scheme I and Scheme II event selection 
schemes, respectively. 



scale choices (about 10%). 

In Figs.4(a,b,c), we present the LO and NLO QCD corrected distributions of the 
transverse momenta of the final Z^-boson, ff'^-boson and leading-jet {p^ \pj^ \pt^) for 
the process pp — )• Z^H^j + X with = /xq at the ^/s = 14 TeV LHC, by adopting 
the inclusive and exclusive event selection schemes, respectively. The corresponding K- 
factors {K{pt) = '^"^^^^ I '^dp!^ ) are also shown in these figures. The analogous plots for 
the pp —7- Z^H^j + X process with ^ = /ig at the -y/s = 7 TeV LHC are presented 
in Figs.5(a,b,c). In these 6 figures, the curves by applying the inclusive event selection 
scheme demonstrate that the NLO QCD corrections generally significantly enhance the 
LO differential cross sections for the pp — t- Z^H^j + X process, especially in the range 
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(c) 



S 3.0 



PP^ Z°H°+jet+X 



Vs=14TeV 



LO 

NLO(I) 
NLO(Il) 
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Figure 4. The LO, NLO QCD corrected distributions of tlic transverse momenta of the final 
particles and the corresponding K-factors {K{pt) = ^^^7^/^^) for the pp Z'^H^j + X process 
at the LHC with jj, = po and y/s = 14 TeV. (a) Z°-boson, (b) Higgs-boson, (c) final leading-jet. 



of ',Pt < 150 GeV, while the curves by adopting the exclusive event selection 

scheme show that the NLO QCD corrections increase the LO differential cross sections 
only in the low transverse momentum ranges of p^. ,Px < 40 GeV and pf, < 70 GeV. 
We also see from Figs.4(a,b,c) and Figs.5(a,b,c) that the NLO QCD correction does not 
obviously change the shape of the LO distribution of the final particles by using both the 
inclusive and exclusive event selection schemes. 
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Figure 5. The LO, NLO QCD corrected distributions of the transverse momenta of the final 
particles and the corresponding K-factors {K{pt) = ^^^7^/^^^) for the pp — )■ Z'^H^j + X process 
at the LHC with /i = /.jq and \/s = 7 TeV . (a) Z^-boson, (b) Higgs-boson, (c) final leading-jet. 



For the further analysis of the uncertainty due to the scale variation, we take another 
scale choice, i.e., by implying the scale //i = \{E^ + Ei^ + Yljet^T^) which is relevant 
to the transverse energies of final particles as the factorization/renormalization scale ^ = 
fif = [46], and compare with the case of /x = //q- The LO and NLO QCD corrected 
transverse momentum distributions of final Z'^-boson, ff'^-boson and leading-jet, and their 
corresponding K-factors with /i = ^1 at the = 14 TeV and ^/s = 7 TeV LHC, are 
plotted in Figs.6(a,b,c) and Figs.7(a,b,c), separately. In comparison with the corresponding 
curves in Figs.4(a,b,c) and Figs.5(a,b,c), we can see the transverse momentum distribution 
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Figure 6. The LO, NLO QCD corrected distributions of the transverse momenta of the final 
particles and the corresponding K-factors {K{j)t) = ^^^7^/^^^) for the pp Z^H^j + X process 
at the LHC with ^ = /ii = ^{E^ + + Ejet-^r*) and ^/i = 14 TeV. (a) Z"-boson, (b) 
Higgs-boson, (c) final leading-jet. 



curves for two scale choice are quite similar, but the curves of Jr(pT)-factor for /i = /ii are 
somewhat more stable than the corresponding ones for /i = /io- It seems that the scale 
choice of the phase space dependent scale // = //i is more appropriate than the fixed one 
/X = /io for the Z^H^ + jet production process at the LHC. 
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Figure 7. The LO, NLO QCD corrected distributions of the transverse momenta of the final 
particles and the corresponding K-factors {K{pt) = ^^^7^/^^^) for the pp — )■ Z'^H^j + X process 
at the LHC with /i = /ii = \{E^ + E;^ + J2jet *) and = 7 TeV. (a) Z"-boson, (b) Higgs- 
boson, (c) final leading-jet. 



4. Summary 

In this work, we present the a full treatment of the next-to-leading order QCD corrections 
to the observables of the pp — ?• Z^H^j + X at the early (^/s = 7 TeV) and future (-y/s = 
14 TeV) LHC. We investigate the dependence of the LO and the NLO QCD corrected 
cross sections on the renormalization/factorization scale, and study the LO and NLO QCD 
corrected distributions of the transverse momenta for the final particles {Z^, H^, leading- 
jet) by adopting the inclusive and exclusive event selection schemes. By taking fi = fiQ = 
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{mz + mH)/2, rriH = 120 GeV and the jet constraint of pf. > pj^j = 50 GeV, we find 
the K-factors for the total cross section have the values of 1.183 (0.988) at the ^/s = 
14 TeV {^/s = 7 TeV) LHC by adopting the inclusive and exclusive schemes, respectively. 
Our numerical results also show that the NLO QCD corrections obviously modify the LO 
integrated and differential cross sections, and significantly reduce the scale uncertainty of 
the LO cross section by adopting either the inclusive or the exclusive schemes. We find 
also that a scale choice like = fii = \{Ej^ + Ej^ + E'ljf*) seems to be more appropriate 
than the fixed scale ^ = fiQ. 
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A Some expressions of the amplitudes at LO and NLO 

The partonic processes of the Z^H^ + jet production at the LHC are related to the 
amplitudes of partonic processes q{pi)q{p2) — ^ {p?,)H^ {VijoiPb) and q{q){pi)g{p2) — ^ 
Z'^{p3)H^{p4)q{q){p5). The summations of the squared amplitude over all spins and colors 
for the q{pi)q{p2) Z'^{p^)H^{pi)g{p^), q{pi)g{p2) Z^iP3)H^{Pi)q{P5) and q{pi)g{p2) 
Z^{p^)H^{p4^)q{p5) partonic processes at the LO are expressed separately as 



-^LO = 2 ft ^ \2 2m u "^^ 2pi ■P5P2-P5- m,pi ■ P2) 



spin 



452e4m2^(32s2,-24s2, + 9) 

+2pi ■ P3{2p3 ■P5P2- P5)] 



^4vS^ml[{p3 + p^y - m\Y{p2 -P3- PiY 
{2pi • Pi[{2p2 ■ Pzf - 2p3 ■PiP2-p?,- in%p2 ■ Pi] + 2pi ■ p4[-m%p2 ■ Pi 
+P2 ■ P4{mjj - 2m| - 2p2 ■ p^)] + Pi ■ P2[-^{P2 ■ Pzf + 4p3 ■PAP2-Pz + 'm\{'m% 
^2^. 4g.Vm^^(32g2^-24sf,, + 9) 

+ 2p3-pi)\] + — g-^ 2— ■ "2 2l27 ^2? ^ 

^c^s'^m^z[{p3 + p^y - m^z\ {.P2 - Pb) {P2 -p-i-PiY 
{m%{pi ■P3 P2-P5-P1- PaP2 ■ Pa) + Pi ■ P5[mzP2 ■ Pa - (2p3 • Pa + rn%)p2 ■ ps 
+2(?'2 • Psf] + P2 ■ P2[in%(2p2 ■ Pa + 4p2 • P5 - 3p3 • P5 - 3p4 • Ps) + 2(p3 • Pa 
-P3 ■ P5 + 2m|)p2 • Ps] - 2pi ■ pa{p2 ■ Psf + 2pi ■ Ps P2 ■ P3 P2 ■ Pa 
-2pi -ps P2- P3 P2-P5 + 2p3 - Pa Pi- P3 P2-P5 + 2p3 ' Pb Pi ' Pa P2 ■ P3 
-2p3-p5 pi-ps P2-Pa}, (A.l) 



E l-^iol' = \Mfo\'{P2 ^ -P5}, (A.2) 

2olor 
spin 

\MIU' = iMfoPfri ^ -P5,P5 ^ -P2,P2 ^ Pi}, (A.3) 



color 
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where the summations are taken over spins and colors of all particles involved. 

The real emission partonic processes for the Z^H^ + jet production at the LHC are 
related to the amplitudes of — )■ qqggZ^H^ and — >• qqq'q'Z^H^ by crossing symmetry. 
The amplitude of the — )• q{pi, ci)q{p2, C2)g{p3, C3)g{p4, C4)Z^{p5)H^{pq) partonic process 
at the tree-level is expressed as below. 



MiO^qqggZ^H'') = 

Tct.T^%,e;^ {p,)s;^ ips)s;.^ {p,)Czzh _ 

(A.4) 



The amplitude of the — ^ q{pi-,ci)q{p2,C2)q'{p3,C'i)q'{p4,,Ci)Z^{p^)H'^{pQ) partonic 
process at the tree-level can be expressed as 



M(0 ^ qq(^q'Z^H^^ 



r<ii\ ( \ I '^ciC2'^C3C4,^*^i{P5)CzZH9i/ii/2 ^^ui I \-( \ 

^fe)C;V.^(-?'4)^(^'l)^^,''<^z(^^ +]/6)C,^|3^(-p2), (A.5) 

where pj and Cj are the momentum and color of the i-th particle, x is the color index of 
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the gluon/quark propagator and 
i e mw 



CzzH 



cwsw 
i e 
cwsw 



Cqqg = -i gsl^ 



M21 



(A.6) 

(A.7) 
(A.8) 



In above equations /| and Q„ are the weak isospin and the charge of quark (7, respectively. 
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